Renormalized Field Theory of Polyelectrolyte Solutions 
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We present a field-theoretic renormalization group (RG) analysis of a dilute solution of flexible, 
screened polyelectrolyte chains (Debye-Hiickel chains) in a polar solvent. We perform a virial expan- 
sion to calculate the scaling dependence of the Osmotic pressure on the chain density, chain size and 
screening length. The formalism is developed for general chain distributions and we give explicit 
forms for a monodisperse chain size distribution. We utilize the fact that a Debye-Hiickel chain 
may be mapped to a local field theory with two length-scales requiring the calculation of scaling 
functions as well as exponents to fully describe its scaling behaviour. 
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Polyelectrolytes are polymers with ionizable monomers 
which in polar solvents, such as water, dissociate into 
charged polyions and small 'counter-ions' of opposite 
sign. They are of widespread importance with many ap- 
plications in physics, biology and chemistry Q. Their 
solublility in polar solvents means that they also have 
many industrial applications (e.g as super-absorbers). 
Typical examples are DNA and sulphonated polystyrene. 

Polyelectrolytes are quite poorly understood, mainly 
because it has been difficult to deal with their long- 
range coulomb interactions theoretically. In addition, 
screening, complex-formation, counter-ion condensation 
and salt and the comparable length-scales of their in- 
teractions make comparison between theory and experi- 
ment difficult. Neutral polymers on the other hand, have 
been extremely succesfully described by scaling ideas and 
renormalization group theories |^| . We recently proposed 
a new approach to the renormalization group applied to 
the study of weakly charged flexible Debye-Hiickel (DH) 
polyelectrolytes in solution || . 

We first review some of the results and the philoso- 
phy of the previous paper. We found that the 'range' 
of the interaction, the Debye screening length, was an 
important new length-scale independent of the system 
size. This is fundamentally different to the systems usu- 
ally studied in conventional statistical mechanics with 
short-ranged interactions where the size of the system or 
volume is the only relevant macroscopic (low q cut-off) 
length-scale and all extensive quantities scale with this 
volume. We obtained non-trivial scaling of quantities 
like the end-to-end distance and radius of gyration of iso- 
lated chains with the number of monomers and with the 
screening length . This second scale also means that 
naive scaling arguments are liable to break down. These 
complications are not new: it has long been understood 
H that taking the thermodynamic limit is problematic 
in systems with long-range Coulomb interactions. 

In the literature regarding flexible chains with un- 
screened coulomb interaction, it is generally agreed that 
the chains are in a 'rigid' regime in d = 3. It is then 



tempting to suggest that even if one has screened inter- 
actions there is a length-scale much below that of the 
screening length where the interaction is effectively un- 
screened and thus the chain would be in a rigid state and 
not amenable to description by a flexible chain model. An 
important point has been overlooked in this analysis: the 
physical arguments leading to the rigid regime depend on 
the region of the chain being considered being 'extensive' 
i.e. of the order of the system size, but if one is dealing 
with screened interactions the region of the chain where 
the interaction is 'effectively' unscreened is by definition 
not extensive. There are fluctuations happening on much 
longer length-scales. In experiments polyelectrolytes are 
found to be extended but not rigid jj) . We conclude that 
there is, in principle, no problem with a flexible model 
being used to describe realistic screened polyelectrolytes 
as long as one is studying the long-wavelength behaviour. 

Another critique of flexible chain models of polyelec- 
trolytes is related to the phenomenon of 'electrostatic 
rigidity'. This is due to the repulsion of regions of the 
chain far from each other. Chains which are much longer 
than the screening-length can be partitioned into flex- 
ible coils or 'Pincus blobs' which align into a rod-like 
blob-pole Such long chains are outside the regime 
of our theory. In fact it has recently been shown that 
the number of monomers of polyelectrolyte required for 
the polyion to be long enough to be in this rigid Pincus 
blob-pole in normal solvents is several orders of magni- 
tude larger than any physically or computationally real- 
isable polyelectrolytes The indication seems from 
experiment and simulations that the viable and inter- 
esting region of polyelectrolyte behaviour is that when 
the screening-length and the chain size are comparable 
and large (macroscopic) - exactly that region being ad- 
dressed by our theory. The first letter addressed single 
chain properties. In this Letter, we apply this formal- 
ism to the calculation of contribution of the polyions to 
the osmotic pressure at non-zero chain concentrations, a 
problem of experimental relevance. 

In the presence of counter-ions and/or salt the coulomb 
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interaction becomes screened. The simplest and most 
widely used model of a screened polyelectrolyte which we 
use as the starting point for our analysis is a solution of 
flexible chains in d dimensions with a Debye-Hiickel (DH) 
or Yukawa interaction between charged monomers 0. 

This is given by V(r) = g J j^y-^fz^ which has the 



familiar form V(r) 



fie 



- K l r l/|r| in d = 3. The chains 



are described by a set of vectors {i"j(s)} parametrised 
by the arc- length s and label j. We define a l Kuhn- 
segmentf length of the chain, £ and Lj — £ 2 Nj is the 
size (area) of the jth chain, where Nj is the number of 
segments on that chain. The Kuhn-length here refers to 
the distance between charged monomers. In d = 3 the 



coupling constant g oc As = (3 



47T££o 



is proportional to 



the Bjerrum length, A^ which measures the strength of 
the bare coulomb interaction, where 1//3 = kgT. In wa- 
ter, \ B = 7.14A. As usual ee is the dielectric constant 
of the solvent, k B is the Boltzmann constant and T the 
temperature. The screening length or range of the inter- 
action is given by l//t and k is a function of the density 
of screening ions and the dielectric properties of the solu- 
tion. The flexible chain DH model is thought to be valid 
only when the polyelectrolyte is weakly charged and in a 
non-zero ion density 

We calculate the contribution of the polyions to the 
osmotic pressure for a solution of chains of N monomers 



and screening length X/k and chain density p p . We find 
the general scaling form 



n p (p p ,iV, k) 
p p k B T 



l + Fi(p p £ d /R d ,NK 1 /"<t>) 
1 + F 2 {p p i d N dv 'i>, R 2 N 2v <»). 



(1) 



This scaling law is the screened polyelectrolyte equivalent 
of the des Cloizeaux scaling law |l0| for neutral polymers. 
From the virial expansion we find explicitly that for a 
monodisperse distribution of chains of N monomers 



3D Pp £ d 



16k 



(2) 



where e = 6 — d, v$ = 2/(d — 2) and D is a non- universal 
constant. P(x) is a complicated function and we have 
defined a dimensionless screening variable measuring the 
screening-length in terms of the Kuhn-length, R = k£. 
These are the main results of our letter. What follows is 
a brief outline of the calculation; details will be presented 
in a longer future publication |llf| . 

We find it convenient to work in the Grand Canonical 
Ensemble of chains with chemical potentials controlling 
the monomer and chain fugacities. The grand canonical 
partition function is given by 



Z[{a}} = 1 



w 1 pOQ 



.dL p e a{Ll) 



,L p ,k] 



(3) 



where <j(L) is the fugacity of a polyion of size (and not necessarily a linear function of) L and 

Z[ Ll , . . . , L p , K ] = J p 2>M.)]e-* ^ *(**/*>'-»E^ C So"' ^vCr.W-^O). 



(4) 



Following de Gennes Q| and des Cloizeaux and performing an Hubbard-Stratonovich transformation p3| , pT| , 
this may be mapped to a local-field theory with two sets of fields with a 'magnetic' field coupled to one set of fields. 
The grand canonical partition function, Z[{er}] is given by 



J 3=1 



(5) 



where the ^ (r) are n — > component (polyion) fields 
with O(n) symmetry and ip(j) is a scalar (coulomb) field. 
The tj are the polyion monomer fugacities, 1/k 2 is the 
screening length andW= (27r) _1 / 2 (det{V 2 + K 2 }) 1 / 2 is a 
normalisation factor. The coupling is given by u = li^Jg. 
The resulting field theory is different to that from neu- 
tral polymer solutions |2j where we have only one type 
of field and the interaction potential is quartic f r (<fr 2 ) 2 . 
Here we have two types of field and a cubic interaction 
f r 4> 2 ip [ fust - The upper critical dimension d c = 6 fol- 
lows from a dimensional analysis of the action in equa- 



tion (^). The 'magnetic' field Hj which without loss 
of generality we align along the first component m of 
the polyion fields, is related to the polyion fugacity, 
e a(L *> = f(L) = Y. i= i e ~ Lt '\ H h li is useful t0 rewrite 
this in the continuum form f(L) = J ■^- i e~ Lt \h 2 (t) mak- 
ing the association ■^rh 2 (t i ) = H 2 . 

We want to calculate the contribution of the polyions 
to the osmotic pressure, n p given by 

\\ = k B T\ogZ[{a}}. (6) 

There is a simple relation between the polyion fugacity 
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<j(L), the monomer density p p and the chain-size distri- 
bution. The number density of of polyions of size L is 
given by 



SH 2 



Sa(L) 'df(L) 

For a monodisperse solution of polyions of size Lq, 
C P (L) = p p S(L-L ). 



(7) 



(8) 
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We can easily generalise our results to arbitrary distri- 
butions by changing the form of C P (L). We perform a 
virial expansion of the osmotic pressure calculating the 
first and second virial coefficients using renormalized per- 
turbation theory. This entails doing a double expansion 
in H and u. A typical term in our perturbation series 
to m-th order H and l-th order in u (with m, I even) is 
given by 



4 (0) . . . 4>) m (0) / &(-k - q) • <MkMq) 



k.q 




H=0 



where C refers to connected diagrams. The series for 
lip can be calculated using standard diagrammatic meth- 
ods and we show the diagrams required to calculate the 
first and second virial coefficients in Fig. 1. 

Renormalization proceeds via the absorption of ultra- 
violet divergences in Z-factors which are sums of poles in 
e = 6 — d. We use a dimensional regularisation scheme 
with minimal subtraction of poles p6j . We can express 
the renormalized parameters in terms of the bare param- 
eters 



u 

K 2 



^/ 2 S^Z u u R 
p 2 (n 2 ) R 

p 2 z^ [tfZ^ + (k 2 ) r z 



= M" 



l Z~ 1/2 H R 
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(9) 



where p 1 is an external length-scale and Sd = Vld/{2ii) d 



with fid the angular part of a d-dimensional volume inte- 
gral. It is important to note the matrix renormalization 
of the polyion fugacities. Because the chains renormal- 
ize independently, the RG flow and fixed point values of 
u R * are the same as for the single chain ||. In particular, 
working at the non-trivial D-H fixed point, we calculate 
the exponent v$ = 2/(d — 2). 

We define the renormalized osmotic pressure 
IL p /k B T = p, d IL R (p R ,L R , (k 2 ) r ,u r ). We can solve it's 
RG equation to give us the general scaling form of Hp. 
We extract the scaling behaviour of the osmotic pressure 
by writing IT^ in terms of flowing variables, fixing the 
RG flow parameters by a matching condition and then 
transforming back to the bare quantities. The result is 
equation (0). The renormalized osmotic pressure can 
also be calculated to 1-loop order in a virial expansion 
(see Fig. 1) as 



n, 
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dt R (h R (t R )) 2 (h R (t R )f 
2ni 2 2 



a 2 (t E 
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0((h R ) 



(10) 



where a\(t, k) is the contribution from the 2-point diagrams and ci2(ti,t2, n) is the contribution from four-point dia- 
grams. Using equations (0), (||) and (JTTj), performing multiple inverse Laplace transforms, we find for a monodisperse 
solution of polyions of size L, 



n 



p P 
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where A\(L,k) 



I 5f? e * ia i(*' K ) and A 2(Li,L 2 ,k) 



Ax *(L,k) = A^iL*}/"*) = 



f rftii*a e tiii e t2i2 a2 ^ l!f2)K ). -yVe find explicitly that 



l + 3ee LK ""^(Lk 1 /^ 



2-iri 2-xi 

and A2(L,L,k) 
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(Lk 1 ^* + \){\n{LK 1 / v i>) +7) - \Lk x I v * - it 



,2Lk 1/v <I> 



Taking into account dimensional factors we obtain 
equation @ with P(x) is the 0(e 2 ) part of A^ 2 (x)A 2 {x). 
It is reassuring that the explicit form satisfies the general 
solution in equation (|I|). 

We compare our result, equations with the sim- 
ulations performed by Stevens and Kremer in d — 3 on 
DH polyelectrolytes [Q. They calculate the scaling be- 
haviour of the polyion contribution to the osmotic pres- 



sure with the monomer density p m — p p N. Their data 
which is for salt-free solutions, i.e. from charge neutral- 
ity k — (47rAsp m ) 1 / 2 , has the osmotic pressure as inde- 
pendent of chain length for both dilute and semi-dilute 
monomer concentrations. For small concentrations they 
find PH P ~ p m whilst for large concentrations they ob- 
serve the neutral polymer scaling (3Ii p ~ Pn( 4 - This 
is consistent with our general scaling law in equation 
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(U) evaluated at d = 3 with F± given by F±(x,y) = 
— 1+ f\xy + fax y + g(x, y) where fi and f% are non- 
universal constants. The /i term (which is the source of 
the small p m behaviour) has exactly the form as that in 
the leading term of our virial expansion, equation (|^). 
The term g(x, y) measures the N dependent cross-over 
between the two regimes. We point out the interesting 
fact that from the structure of the field theory, there are 
three fixed points governed by the ratio c = £ 2 k 2 N ||. 
Because n ~ pin, 2 , by varying p m , we cross-over from the 
polyelectrolyte (c ~ 1) to the SAW fixed point (c >> 1). 
Our calculation explicitly captures this cross-over. The 
third fixed point (rod) occurs only at p m = 0. One may 
wonder why although we have performed an expansion 
in e = 6 — d and e — 3 in 3 dimensions, we get rather 
good agreement with the simulation data. We suggest 
that this is due to the symmetry identity of the vertex 
functions of the field theory which leads to the exact cal- 
culation of the exponent v$ |Q . Our scaling law can also 
be used for polyelectrolytes in a salt solution but there 
the leading behaviour of k is not related to the monomer 
density Q. 

It is also interesting to analyse the form of the second 
virial coefficient from equation (|^). In the theoretical 
work on neutral polymers, this quantity is used to de- 
fine the 'overlap' variable which usually has the form (in 
d = 3) (i? 2 ) 3 / 2 ~ N 3ui where (R 2 ) is the mean end- 
to-end distance and vi — 0.588 is the size scaling expo- 
nent of the Edwards SAW. Here we find that the overlap 
variable is rather complicated function of the screening 
length, 1/k and the chain size N. This is reasonable for 
a system with long-range interactions as the 'effective' 
size of the particles strongly depends on the range of the 
interaction as well as the physical size of the object. We 
see that the overlap variable can in principle be some- 
what larger than the chain size. This is also consistent 
with the observations in simulations |l7[ . 

In summary, we have calculated the scaling behaviour 
of the contribution of the polyions to the osmotic pres- 
sure using field theoretic methods. We obtained a gen- 
eral scaling law of the osmotic pressure. To our knowl- 
edge, this is the first systematic attempt to calculate the 
osmotic pressure of a solution of Debye-Hiickel chains 
without uncontrolled approximations ||. We also found 
that we can explain all the main features of the osmotic 
pressure measured in some recent simulations of screened 
polyelectrolytes. At a more fundamental level, this calcu- 
lation explicitly shows the subtleties of performing field 
theory with two relevant and independent length-scales 
and provides another interesting application of field the- 
oretic methods to soft condensed matter systems. 
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Stevens and T. Vilgis. We particularly thank K. Kremer 



and T. Vilgis for a critical reading of the manuscript. 
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FIG. 1. The diagrams for (a) di(t, k) and (b) 02(ti,i2, «)■ 
A straight line corresponds to a free polyion propagator 
(q 2 + tj) _1 , a wavy line a coulomb propagator (q 2 + k 2 ) -1 
and a vertex to u/2. We integrate momentum over all closed 
loops 
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